STUDIES OF A 2-CHAIN SPIN LADDER 
WITH FRUSTRATING SECOND NEIGHBOR 
INTERACTIONS 

Zheng WeihongQ, V. Kotovi, J. OitmaJ 
School of Physics, The University of New South Wales, Sydney, NSW 2052, Australia. 

(February 1, 2008) 



Abstract 

The Heisenberg model on a 2-chain spin-^ ladder with frustrating second 
neighbor interactions is studied by using series expansions about the Ising 
and dimer limits, numerical diagonalization, and by Abelian bosonization 
analysis. The phase diagram is determined, and pair correlations and the 
complete dispersion relations for the triplet spin-wave excitations are also 
computed. 

PACS Indices: 75.10.-b., 75.10J., 75.40.Gb 
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I. INTRODUCTION 



Heisenberg spin ladders have been the subject of intense theoretical and experimental 
research in recent years. It is by now well established that single chain Heisenberg anti- 
ferromagnets with integer spin have a gap in the excitation spectrum, whereas those with 
half-integer spin have gapless excitations. The former have a finite correlation length, while 
for the latter it is infinite with the spin-spin correlation function decaying as a power law. 
For S = ^ Heisenberg spin ladders , the systems with an even number of legs have an 
energy gap, short range correlations and a "spin liquid" ground state. On the other hand, 
the systems with odd number of legs have gapless excitations, quasi long range order, and 
a power-law falloff of spin-spin correlations, similar to single chains. 

In this paper, we study the Heisenberg model on a 2-chain ladder with second neighbor 
interactions, via Ising and dimer expansions, diagonalization of finite systems, and Abelian 
bosonization. All of the work is at T = 0. The motivation for studying such a system 
is twofold. Firstly we wish to explore the effect of frustration on the properties of 2-leg 
s = 1 ladders. Weak frustration is not expected to change the physics of the gapped system 
qualitatively. However strong frustration may change the nature of the ground state at 
some critical coupling ratio, corresponding to a T = phase transition. Furthermore in 
other systems, such as the Ji — J2 chain, frustration itself leads to the creation of a gapped 
dimer phase, and thus in the present case there is the possibility of observing the competition 
between two independent gap yielding perturbations. The second reason for studying such 
a generalized ladder system is the possibility that in real 2-leg ladder materials significant 
second neighbour interactions will be present. 

We write the Hamiltonian of our system as 

i 

+■^2(81,4 ■ 82,4+1 + 82,4 ■ Si,j+i)] (1) 

where S;,j denotes the S = 1/2 spin at site i of the Ith chain. J^^ is the interaction between 
nearest neighbor spins along the chain, J± is the interactions between nearest neighbor spins 
along the rungs, J2 is the interactions between the second neighbor spins. This is shown 
in Figure 1(a). We denote the ratio of couplings as yi = Jj_/J^^ and 1/2 = J2/J//- In the 
present paper all couplings are taken to be antiferromagnetic (that is, J^^, J±, J2 > 0). 

The system has an interesting symmetry property: If one exchanges the couplings J^^ 
and J2, one can recover the original Hamiltonian by exchanging two spins along the rungs 
at even sites, that is, the Hamiltonian will be invariant under exchanging the couplings J^^ 
and J2: 

HiJ^^,J^,J2) = H{J2,J^,J^^) (2) 

Therefore we only need to study the case of 1/2 < 1, and the system with j/2 > 1 can be 
mapped into system with ?/2 < 1 through the identity in Eq. (Q). 

There has been some previous work on this system. Of course for J2 = we recover 
the usual 2-leg ladders |T|-§[, which has gapless excitations only when J± = 0. The case 
1/2 = 1 has special properties. Bose and Gayen first pointed out that this system has 
an exact dimer state: a state in which every pair of spins along the rungs form a singlet, 
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and this perfect dimer state is the ground state for large enough values of yi. Xian 
performed a systematic study of this system via a microscopic approach based on a proper 
set of composite operators, and found that the Hamiltonian consists of two commuting parts: 

H = ^2 ^-L^l,« ■ ^2,i + J//Pi ■ Pj+1 (3) 
i 

where Pj represent the effective spin-1 operator at site i. The first part, J_LSi^j ■ 82,,, is 
trivial, with no interaction between different rungs and with a gap J± between the singlet 
and triplet states. The second part, J //Pi ■ Pj+i, is similar to the spin-1 Heisenberg chain 
with each rung in the ladder corresponding to each site in the spin-1 chain. Because the two 
part of H commute it follows that there is an eigenstate of H with eigenvalue: 

Eg = (J//eo + JjA)N/2 (4) 

where Cq = —1.40148403897(4) is the ground-state energy per site of the spin-1 Heisenberg 
chain . The eigenvalue for the state with singlets on all of the rungs is 

Eg = {-3JjA)N/2 (5) 

Therefore one can easily determine the transition point i/c = 1.40148403897(4). At this 
crossing point the singlet-singlet gap is zero. Kitatani and Oguchi ^ had independently 
considered the existence of a state of singlet dimers for this system for 1/2 = 1, and had 
obtained the same results as above. 

This paper studies the general case where exact results are not known. In Section 2 
we report results obtained by Ising and dimer expansions at T = 0. We compute the 
ground state energy and the complete spin-wave excitation spectra, as well as obtaining more 
accurate estimates for various quantities for the simple case J2 = 0. We also carry out finite 
lattice diagonalization studies for systems of 16 and 24 spins, paying particular attention 
to the behaviour of spin-spin correlations. Section 3 describes an analytic approach, based 
on Abelian bosonization. When comparisons are possible the various methods all provide a 
consistent picture. Finally in Section 4 we provide a summary and draw some conclusions. 

II. NUMERICAL RESULTS 

We report here the results of numerical studies, using both series expansions and diago- 
nalizations. These methods provide results over the whole phase diagram, and give the most 
detailed picture of the behaviour of this system for general values of the couplings. 

Our previous series work on the ladder system used expansions about both the Ising 
limit and about a fully dimerized state, and we refer to that paper and references therein 
for technical details. 

A. Ising expansions 

To construct an expansion about the Ising limit for this system, one has to introduce an 
anisotropy parameter x, and write the Hamiltonian as: 
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H = Ho + xV (6) 

where 

= ('^M'^M+i '^2,i'S'2,i+i) + J:LSliS2^^ + J2{,Sl^^Sl,ij^i + 5'£jS'^^j_|_;^)] (7a) 

i 

T/ _ Y^r J ( ox QX j_ oy oV 1 ca; cas- j_ oy oV \ \ J ( \ qv Qy \ 

+ J2(S'f^jS'2_j+i + + 'S'2iS'J^^j_,_i + 5'|_iS'f^j_,_;^)] (7b) 

The hmits x = and x = 1 correspond respectively to the Ising and isotropic Heisenberg 
hmit, the latter being the case of primary interest. 

Since Hq is taken as the unperturbed Hamiltonian we need to identify the ground states. 
There are three of these, shown in Figure 1(b). We refer to these as the Neel state (N), the 
ferromagnetic chain state (F) and the ferromagnetic rung state (R). Their regions of stability 
are indicated in Figure 2. The line 1/2 = 1, > 2 is a boundary between N and F states. 
The operator V is treated as a perturbation: it flips pairs of spins on neighbouring sites. The 
quantum fluctuations represented by V will of course result in much more complex ground 
states. Since the system is effectively 1-dimensional no true long range order can exist, even 
at T = 0. 

As in our earlier paper to overcome a possible singularity at x < 1, and to get a 
better convergent series in the Heisenberg limit, we add to the Hamiltonian the following 
staggered fleld term 

AH = til-x)J2i-iy^'Sl, (8) 

i,i 

for the expansion about the Neel state, or the following fleld term 

AH = t{l-x)J2i-iySl, (9) 

for the expansion about the ferromagnetic rung state. AH vanishes at the isotropic limit 
X = 1. We adjust the coefficient t to get the most smooth terms in the series, with a typical 
value being t = 2. 

The series expansion method has been previously described in several articles |]8HlO| , and 



will not be repeated here. We have developed Ising expansions about both the Neel and 
ferromagnetic rung states for the ground state energy per site Eq/N and the triplet spin- 
wave excitation spectrum for different ratios of couplings yi and 1/2 and (simultaneously) for 
several values of t up to order x^^. The resulting series are available on request. These series 
have been analyzed by using integrated flrst-order inhomogeneous differential approximants 



and Fade approximants |11]. We will discuss these results later in this section. 



B. Dimer expansions 

In the limit that the exchange coupling along the rungs J± is much larger than the 
couplings J . and J2, the rungs interact only weakly with each other, and the dominant 
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configuration in the ground state is the product state with the spins on each rung forming 
a spin singlet. The Hamiltonian in Eq. (|l|) can then be rewritten as, 

H/J^ = Ho + {l/yi)V (10) 

where 

i 

(11) 

V = X/[(^l.« ' 8l,i+l + ^2,i ■ 82,1+1) + 1/2(81,1 • 82,4+1 + 82,4 • 8i,i+i)] 

i 

We can obtain an expansion in (l/yi) by treating the operator Hq as the unperturbed 
Hamiltonian, and the operator V perturbation. 

We have carried out the dimer expansions for the ground state energy up to order (l/?/i)^^ 
and for the triplet excitation spectrum up to order (l/yi)^'^ for different values of 2/2 • The 
resulting series for some particular values of y2 are listed in Table I. The rest of the series are 
available on request. In our previous paper 0], the series for the case 1/2 = were computed 
up to order for the ground state energy and up to order (l/yi)^ for the triplet 

excitation spectrum. Our present results agree with and considerably extend these previous 
results. Again, we use integrated first-order inhomogeneous differential approximants and 



Pade approximants [11] to extrapolate the series. 

The ground state energy per site Eq/N is shown in Figure 3, where curves for y2 = 
0,0.2,0.4,0.6,0.8,1 as shown as functions of yi = J±/J^^. The curves (connecting the full 
point symbols) emanating from the left side (small yi) are obtained from Ising expansion 
about the R state. Those emanating from the right side are from the dimer expansion. The 
estimates from the Ising expansion about the N state agree very well with these and are not 
shown separately (except for the case of ?/2 = which are shown by open point symbols with 
error bars). These curves, for any given 7/2, cross at a transition point yic. This corresponds 
to a first order ground state phase transition, resulting from a level crossing. The numerical 
estimate for ?/2 = 1 of yic = 1.40 is in good agreement with the exact result of Xian 
and Kitatani and Oguchi discussed above. The locus of the transition points is shown 
in Figure 2, and represents a "gapless line" where the gap between the two lowest energies, 
both of which are singlets, vanishes. Everywhere else the ground state is a non-degenerate 
singlet . However this is not a gapless line in the usual sense of a vanishing gap for elementary 
excitations. 

There are two branches of triplet spin excitations which we denote eg{k), eu{k) corre- 
sponding to symmetric and antisymmetric states with respect to interchange of the two 
chains, k is the wavenumber along the chain direction. On the left side of the gapless line 
in Fig. 2, which can be called a "Haldane" phase, these two branches are independent and 
each has the appearance of a simple cosine dispersion curve, symmetric about k = 7t/2 and 
having a finite gap at the minimum energy points. However on the right side of the gapless 
line, the "dimerized region", the branches are no longer independent, being related by 

eg{k) = 6^(71 - k) (12) 

and having a complex form, without the symmetry about k = tt/2 found in the Haldane 
region. This relation can be understood as follows. Consider an initial down spin at site j 
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flipped to create an excitation. To form a Bloch state tliis excitation will couple to other 
downs spins. In the Neel or small ?/2 region this involves even sites on the upper chain and 
odd sites on the lower chain leading to a shift of vr in the wavenumber. The dispersion curves 
are shown in Figs. 4-5 for yi = 1 for particular values of ?/2- The triplet gap appears to be 
nonzero throughout the phase diagram, except at yi = y2 = 0. 

Since the dimer expansion carried out here is much longer than our previous calculations 
1^, we can make more accurate estimates for the normal 2- leg ladder without frustration 
(J2 = 0). The ground-state energy is estimated to be 

Eo/{NJ^^) = -0.578043(2) at J^^ = . (13) 

In Fig. 6 we show the dependence of the triplet excitation gap eu(7r) on J±/{J± + J^^), in 
particular, they are estimated to be 

e„(7r)/J^^= 0.5028(8) at J^^ = ; (14) 
e«(vr)/J±= 0.405(15) at J > . (15) 

These results agree very well with the recent Quantum Monte Carlo (QMC) results of 



Frischmuth et al. [12] and Greven, Birgeneau, and Wiese 



C. Exact diagonalizations 

In order to obtain a more complete picture of the energy spectrum, the correct assignment 
of spin quantum number S to different levels, and the variation of spin-spin correlations 
throughout the phase diagram we have carried out exact Lanczos diagonalizations for ladders 
with = 16 and = 24 spins. The finite lattice corrections appear to be quite small and 
the results are believed to be representative of the thermodynamic limit. 

In Fig. 7 we show the energies of low lying energies for N = 24, for a scan through the 
phase diagram (Fig. 2) at 7/2 = 0.8. The solid lines represent singlets (5* = 0), while triplet 
levels (5* = 1) are shown as dashed lines. The existence of the transition point discussed 
above, the finite singlet-triplet gap, and other level crossings can be seen. There is a region 
around the transition point where the lowest triplet lies above the lowest two singlets. The 
location of the transition point, while less accurately resolvable, is totally consistent with 
the series results given above. Fig. 8 we show some pair correlations, for the same scan 
through the phase diagram at ?/2 = 0.8. The correlations show large discontinuities at the 
transition point yi ~ 1.24. The sign of various correlations is consistent with ferromagnetic 
rung (R) type order for yi < 1.24 and Neel type order for yi > 1.24, in accordance with the 
classical ground states in Figure 1. We note also that further neighbour correlations become 
very small for increasing yi, consistent with a dimerized phase. 



III. WEAK COUPLING ANALYSIS 

In this section we present analysis, based on Abelian bosonization, which can determine 
the phase boundary between the Ferromagnetic rung (Haldane-type) and Neel ground states. 
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Our considerations are valid in the limit of small J±_ and J2. Since the procedure is well 
described in the literature 



13-15 



we will give only the basic steps here. 
The spin operators for each chain are transformed, by using the Jordan- Wigner transfor- 
mation, into a system of spinless fermions a„(chain 1) and 6„(chain 2) at half filling. Next, 
since we are interested in the low-energy properties of the model, we pass to a continuum 



description, which was developed for a single chain by Luther and Peschel |T^. The spec- 
trum of the Jordan- Wigner fermions is linearized in the vicinity of the two Fermi points 
±kF = ±tt/2: 



ya[e^^^"^iij(x = na) + e-'^^^^/^ii (x = na)] 



(16) 



where a is the lattice spacing, and ipm, ipiL are slowly varying on the scale of the lattice. 
An analogous transformation is applied to the second chain (with corresponding fields ijj2R, 
ip2L)- In order to write the ladder Hamiltonian (1) in bosonic form, we make use of the 
Abelian bosonization rules: 



1R,1L 



exp[±iiV47r0ii{_iL(a;) 



where a'^ 



~ a^^ is a large momentum cutoff, and: 

. . 1 



hRAL[X) 



[Mx)±e,{x)]. 



(17) 



(18) 



Here 6*1 (x) is defined as the field, dual to 0i(x), i.e. dx9i{x) = ni(x), and ni(x) is the 
momentum field, conjugate to 0i(x). Similar equations describe the second chain. Using 
the above formulas the interchain interactions can be bosonized, with the result: 



H 



V — 

2 



dx 



(19) 



dx 



(7icos(v^87r0_|_) + 5f2COs(v^87r</)_) + 25f3COs(v^27r6'_ 



(20) 



In the above equations, the symmetric and antisymmetric combinations of the fields are 
introduced via v^0± = 0i ± 02? and similarly for the dual fields. The (bare) values of 
the coupling constants in (pOD are Qi = yi — 2y2,i = 1,2,3, and C is a cutoff independent 
constant. The Luttinger liquid parameters in (O) depend on the bare couplings: 



v± = 1 -\ — ± h higher order. 



(21) 



1 — 2— =p — — + higher order. 
IT 2tt 



(22) 



Without the j/2 terms our equations are similar to the ones obtained by Strong and Millis for 
the simple ladder |T^. In order to emphasize that the second terms in (pT]) and ( p2D come 



from the Ising interactions in the two chains, we have introduced the anisotropy parameter 
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S (i.e. JJ^ — > J^^5, and we should set 5 = 1 at the isotropic point). In the hmit yi,y2 — ^ 
the exact form of the above functions are known from the Bethe ansatz solution of the chain 
problem, and, in particular for 5 = 1 we have v± = 7r/2, K± = 1/2 [Q. Away from this 
exactly solvable limit, the expressions (0) and (^21) should be viewed as valid only to lowest 
order in 6, yi and y2, since all lattice renormalization effects have been neglected in passing 
to the continuum limit. We have displayed in (|20|) only the relevant (in renormalization 
group (RG) sense |T^) operators and have neglected all potentially irrelevant and marginal 
ones. The latter contain terms that mix the symmetric and antisymmetric sectors, as well 
as combinations of field derivatives and cosines. 

The scaling dimensions of the three cosine operators in (pO]) are 2i^+, 2K^ and {2K^)~^ 
(corresponding to gi,g2 and g^, respectively). A cosine operators is relevant if its scaling 
dimension is less than two. Thus, in the limit yi = y2 = 0, all operators have dimension one 
and are relevant. For non-zero values of the interchain couplings, one can easily see that the 
most relevant operators are gi and g^. These two couplings thus flow to infinity which signals 
formation of a gap. This strong coupling regime corresponds to a non-zero expectation value 
of < 5*1. 5*2 > [T3|. Whether a ferromagnetic rung (Haldane) or antiferromagnetic ladder 
(Neel) state is realised, depends on the sign of gi and g^. The equations, governing the 
RG flow for gi and g^ are the usual Kosterlitz-Thouless equations [0. Thus, if initially 
(?i(0), 5^3(0) > 0, then gi{l),gz{l) — >■ 00, where / is the RG iteration parameter. This is the 
Neel state, characterized by a gap to triplet excitations. In the opposite limit (71(0), (^3(0) < 
the couplings flow to minus infinity, which is interpreted as a ladder with an effective 
ferromagnetic interchain coupling. Thus we conclude that the transition line between the 
two states is given by the equation gi = gz = yi — 2?/2 = 0. 

Let us note that as yi and y2 increase one should include higher order terms in the 
operator product expansion, which leads to (pOl). The additional terms are still less relevant 
than the ones in (pO]) but they typically couple the symmetric and antisymmetric sectors. 
Thus, the renormalization of gi affects g^ and vice versa. The location of the transition 
line, however, is not affected by this coupling. On the other hand we expect that lattice 
renormalization effects, i.e. inclusion of higher order terms in the lattice spacing, might lead 
to a change of the shape of the transition line, as suggested by numerical simulations (see 
Fig.2.). 



IV. CONCLUSIONS 

We have studied the two-chain antiferromagnetic spin ladder with frustrating second 
neighbour interactions, using a variety of numerical and analytic methods. When the in- 
terchain nearest neighbour coupling is dominant the system is in a gapped "dimerized 
phase" , whereas for dominant second neighbour coupling J2 the system is in a gapped "Hal- 
dane" phase, which can be mapped onto an 5 = 1 chain. These phases have the same 
physical origin within the low energy field theoretic framework but are distinguished by dif- 
ferent behaviour of correlations and are separated by a first order transition line where there 
is a vanishing singlet-singlet gap. The system has a trivial "valence-bond-solid" ground state 
along a line in the phase diagram which separates two types of classical Ising ground state. 
A symmetry of the Hamiltonian allows these states to be mapped onto each other. 
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Using series expansions, about both Ising and dimerized unperturbed states, we have 
computed the ground state energy and dispersion curves for singlet-triplet excitations. The 
latter show a qualitative change in form as the second neighbour interaction changes. Our 
series in dimer expansions are substantially longer than in our previous study for J2 = 0, 
and we obtain very accurate estimates for the ground state and excitation gap, which agree 
very well with recent Quantum Monte Carlo results. We have also computed ground state 
energies and correlations using exact diagonalizations for = 16, 24. These give a further 
physical of the nature of the ground state in different regions of phase diagram. 

Finally we present an analytical study, valid for small and J2, using the technique of 
Abelian bosonization. This gives a picture consistent with the numerical work. 
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FIGURES 



FIG. 1. (a) Three different couplings considered: J n (the dotted lines), J_l (the dashed lines), 
and J 2 (the solid lines), (b) Three different spin orders for the system at the Ising limit: the Neel 
state (N), the ferromagnetic chain state (F), and the ferromagnetic rung state (R). 

FIG. 2. Phase diagram in the (yi, 2/2) plane. The dashed lines show the phase boundary for the 
three different spin orders of the system at the Ising limit. The solid line shows the phase boundary 

between dimerized phase (right side) and Haldane-type phase (left side) for the isotropic system. 
Along this line there is a vanishing singlet-singlet gap. The "x" line (2/2 = 1) 2/i ^ 1.401484) is the 
location of the exact dimer ground states. 

FIG. 3. The ground-state energy per site Eq/ {^N J as a function of yi = J^/J^^ for 
2/2 = 0,0.2,0.4,0.6,0.8,1. The lines in the large yi region are the extrapolations of integrated 
differential approximants to the dimer series, the full point symbols connected by a line are esti- 
mates from Ising expansions about the ferromagnetic rung state, and the open point symbols (for 
2/2 = only) are estimates from Ising expansions about the Neel state. The position of crossing 
indicate transition point with vanishing singlet-singlet gap. 

FIG. 4. The dispersions eu{k) of the spin-triplet excitated states of the 2-chain ladder with 
interchain coupling 2/1 = 1 and 2/2 = 0, 0.2, 0.4, 0.6, 0.7, 0.8, 1. The results for 2/2 < 0.6 are estimated 
from the dimer expansions, and results for 2/2 > 0.6 are estimated from the Ising expansions about 
the ferromagnetic rung order. 

FIG. 5. The dispersions eg{k) of the spin-triplet excitated states of the 2-chain ladder with 
interchain coupling 2/1 = 1 and 2/2 = l, 0.8, 0.7 (shown in the figure from the top to the bottom 
respectively). The results are estimated from the Ising expansion about the ferromagnetic rung 
order. 

FIG. 6. The minimum triplet energy gap e„(7r)/J^^ for the system with J2 = as a function of 
J±/{J±+J//)- The results are estimated from the several different integrated different approximants 
to the dimer series. 

FIG. 7. Low lying energies for a system of iV = 24 spins, for fixed 2/2 = 0.8. Solid (dashed) 
lines represent singlet (triplet) levels respectively. 

FIG. 8. Pair correlations C„ = 4:{SqS^^) in the ground state for = 24 at fixed 2/2 = 0.8. Note 
the discontinuities in correlations at the transitions points 2/ic — 1.24 (which decreases slightly as 
the lattice size N increases). 



11 



TABLES 



TABLE I. Series coefficients for tfie dimer expansion of the ground-state energy per site 



Eo/{NJ±), and the energy gap eu{TT)/J±- Coefficients of " are listed for 2/2 = 0,0.2,0.4,0.6,0.8. 



n 


2/2=0 




2/2 = 0.2 




2/2 = 0.4 




2/2 = 0.6 




2/2 = 0.8 










Ground-state energy Eo/(NJj 


.) 











-3.750000000x10- 


1 


-3.750000000x10- 


1 


-3.750000000x10" 


-1 


-3.750000000x10" 


-1 


-3.750000000x10" 


-1 


1 


0.000000000 




0.000000000 




0.000000000 




0.000000000 




0.000000000 




2 


-1.875000000x10- 


1 


-1.200000000x10- 


-1 


-6.750000000x10- 


-2 


-3.000000000x10" 


-2 


-7.500000000x10- 


-3 


3 


-9.375000000x10- 


2 


-7.200000000x10- 


2 


-4.725000000x10- 


-2 


-2.400000000x10" 


-2 


-6.750000000x10- 


-3 


4 


1.171875000x10- 


2 


-1.920000000x10- 


-2 


-2.548125000x10- 


-2 


-1.770000000x10" 


-2 


-5.981250000x10- 


-3 


5 


8.789062500x10- 


2 


2.880000000X 10- 


2 


— 2.953125000X 10- 


-3 


-1.080000000x10" 


-2 


-5.146875000x10- 


-3 


6 


7.763671875x10- 


2 


4.771200000X 10- 


-2 


1.317346875X 10- 


-2 


-4.452000000x10- 


-3 


-4.292531250x10" 


-3 


7 


-2.682495117x10- 


2 


2.536560000x10- 


-2 


1.777094648x10- 


-2 


3.552000000x10- 


-4 


-3.463308984x10" 


-3 


8 


-1.381530762x10" 


- 1 


-2.345664000x10- 


-2 


1.045650727x10- 


-2 


3.014910000x10- 


-3 


-2.699781797x10- 
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■2 


-4.401953880x10- 
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-3 
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-2 
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-2 
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